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ABSTRACT 

CS Cha is a binary system surrounded by a circumbinary disk. We construct 
a model for the inner disk regions and compare the resulting synthetic SED with 
IRS spectra of CS Cha taken at two different epochs. For our model we adopt a 
non-axisymmetric mass distribution from results of published numerical simula- 
tions of the interaction between a circumbinary disk and a binary system, where 
each star is surrounded by a disk. In particular, we approximate the streams 
of mass from which the inner circumstellar disks accrete from the circumbinary 
disk. This structure is due to the gravitational interaction of the stars with the 
disk, in which an array of disks and streams are formed in an inner hole. We 
calculate the temperature distribution of the optically thin dust in these inner 
regions considering the variable impinging radiation from both stars and use the 
observations to estimate the mass variations in the streams. We find that the 
SEDs for both epochs can be explained with emission from an optically thick 
inner edge of the circumbinary disk and from the optically thin streams that 
connect the circumbinary disk with the two smaller circumstellar disks. To the 
best of our knowledge, this is the first time that the emission from the optically 
thin material in the hole, suggested by the theory, is tested against observations 
of a binary system. 
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Introduction 



Studying the behavior of variability in T Tauri stars surrounded by disks can lead to a 
better understanding of their host stars as well as their disk structure and ev olution. With 



observations obtained by the Spitzer Space Telescope (I Werner et al. 1 120041 ) new detailed 



studies of mid-IR variability have become possible. Several mechanisms have been proposed 
to explain the dif f erent behavior observed in several young circumstellar disks. For instance, 
Muzerolle et al. I J2OO9I ) found that the mid-IR SED of LRLL 31 in IC 348 varied on the scale 
of days, with a characteristic "seesaw" shape, i.e., as the emission decreased at w avelengths 
shortwards of ~ 8.5/im, it increased at longer wavelengths. iMuzerolle et al. I ( 120091 ) proposed 
that this type of variability reflected changes in the height of the inner disk edge. 

Some infrared variability has been explained by c hanging the location of the inner disk 
edge. IR variability was found bv Isitko et al~l (EES) in HD 31648 and HD 163296 in the 
NIR range, (specifically between 1 and 5/zm ), where the changes could be attributed to 
variations in the location of the dust sublimation zone. Modeling under this assumption 
shows that the sublimation wall (boundary between the dusty and the dust-free regions of 
the disk) varies between 0.29AU and 0.35AU for HD 163296. However, any variation of the 
material in the inner region, not limited to changing the location of the inner disk edge, 
could be responsible for the observed NIR variability. 

Some of the intrinsic variability in circumbinary disks comes from the fact that the 
position of the stars is constantly changing along their orbits. Thus, the stars' irradiation 
of the disk also changes, resu l ting i n variation in the observed disk emission. This behavior 
can be seen in iNagel et al. I (120101 ) , where the p redicted variat i on of the dust emission is 



presented for the binary system CoKu Tau/4. IJensen et al. I (120071 ) also show that the 



emission of the binary UZ Tau E is clearly periodical. They explain this behavior with 
emission from material that flows from the circumbinary disk towards the stars. The orbital 
motion of a binary system and its interaction with the circumbina ry disk has been propos ed 
to be responsible for variability on the scale of years for UY Aur (IBerdnikov et al. 1120101 ). 



Espaillat et al. I (1201 ll ) have presented the largest modeling study of mid-IR variability 



in disks around TTS. Most of their sample, composed by transitional and pre-transitional 
disks, show seesaw variability, that was explained as a result of variations in the height of the 
inner disk edge. However, CS Cha, the only known circumbinary disk in the ir sample, shows 



variat ions only in the 10 /im silicate emission band. Based on SED modeling, lEspaillat et al. 
(120111 ) inferred a hole size in the disk of about ~ 38 AU, and this large size implies that the 



inner disk edge is too far from the stars to have a temperature high enough to contribute 
at the ~ 10 /iin emission. In their model, this 10 /zm emission was produced by uniformly 
distributed optically thin dust inside a radius of ~ 1 AU, centered around a single star, and 
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the observed variability was interpreted as due to variations in the mass of this optically 
thin dust. 

CS Cha is an interesting laboratory to compare the observations to simulations of the 
interaction between the circumbinary disk and the binary system. In this paper we construct 
a model that takes in to account results from published numerical simulations of the circ umbi- 
nary disk dynamics ( iGiinther fc Kiev 1 120021 . 12004] ; lArtymowicz &: Lubow I Il994j . 119961 ) , and 
aim to constrain the se simulations by comp aring with the variability observed in the IRS 
spectra presented by lEspaillat et al. I (120 111 ). In particular, we take into account that the 
predicted mass distribution within the inner hole is non-axisymmetric; it is characterized by 
streams connecting the edge of the circumbinary disk with the smaller circumstellar disks 
that surround each star. We find that these streams and the circumstellar disks are the 
main contributors to the emission at 10 /im. Since this material has a very low density, 
we calculate the temperature as a function of grain size assuming each grain is in radiative 
equilibrium with the radiation fields of both stars. 

The paper is organiz ed as follows. § [2] summarizes the observations, already described 



in 



Espaillat et al. I ( 1201 if ). The model is described in section [21 with particular detail given 
to how we estimate the spatial distribution of dust and calculate its temperature. The 
model input parameters, both those taken to be fixed and those adopted as free parameters, 
are described in § HI § [5] describes the results and, finally, § [6] presents the summary and 
conclusions of the paper. 



Observations 



CS Cha was observed at 3 different times. The first spectrum was obtained through the 
IR S Guaranteed- Time Obse r vations (GTO), on J u ly 11th 2005, a n d was previo usly presented 



in 



Espaillat et al. I (120070 1. iKim et al. I (12009( 1. iFurlan et al. I (12009( 1. and iManoi et al 



(120 111 ). The other two spectra were obtained in GO 5 0403 Spitzer prog r am (P hEspaillat) on 
June 1th and 8th in 2008 and were first presented in lEspaillat et al. I (120111 ). The two GO 
observations do not show large differences between them. Therefore, in this work we only 
present the IRS spectrum taken on June 1, 2008. 



Luhman I (120041 ) found an spectral type of K6 for CSC h a; at t hat time, the binarity 



was not established. Using this information, lEspaillat et al. I ( 12007bl ) found a value for the 
extin ction, Ay = 0.8: match i ng th e V-I color of the star to that of a main sequence standard 



from iKenyon fc Hartmann I (11995 ). They take this value to deredden the Spitzer spectrum 
using the Mathis reddening law (IMathis Ill990l ). We adopt the same SED in this work. A 
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distance of d = lQOpc is adopted for this system (IWhittet et al. 1 119971 ) . 



Guenther et al. I (120071 ) found CS Cha to be a binary. They did not report the precise 



orbital parameters, but we use their results as reference. They found a fit with an orbital 



period of 2482 days with a companion of a mass larger or equal than O.1M . lEspaillat et al. 



(120111 ) take a star of 0.91M Q in their single star model. Using the previous restriction on the 
secondary mass and the value of the primary mass in the single star model, we consider two 
cases for the binary system: Case 1 is a set of stars with masses of 0.8M & and 0. 3M M , and 



Case 2 corresponds to 0.9M Q and O.1M stars. For an estimated age of 2Myr (ILuhman 



2004 ) for th e syste m and the masses of the stars, we use the pre-main sequence isochrones of 



Siess et al. I (120001 ) to find the radius R* and the surface temperature T* for each star. Our 



choice is justified by fitting the observed flux at low wavelengths, where the contribution of 
the disk is very small. We decided not to increase the secondary mass value furthe r, because 



i n tha t case, the spectral type of the primary would be far from the one derived by ILuhman 



(120041 ) . In order to argue this, we assume tha t the primary spectral type is close to the 
estimate which takes into acc ount just one sta r (ILuhman 1120041 ) . A comparison between the 
SED of the single star used in lEspaillat et al. I (120111) and the addi t ion of b oth star's spectra 
for Case 1 and 2 is shown in Figure [TJ Finally, lEspaillat et al. I (l2007b| ) estimate a mass 
accretion rate of 1.2 x lO~ 8 M Q yr~ 1 using the UV excess on one star, which we distribute 
between both stars (see § 14.11) . The parameters for both combinations of stars (Case 1 and 
2) are presented in Table [TJ 



The model 



The basic model for CS Cha is similar to the one described in iNagel et al. I (120101 ) . In 
summary, it includes the emission of the inner wall of a circumbinary disk illuminated by 
the two stars of the binary system. It also includes the em i ssion o f optic ally thin dust in 



the inner "hole" of the disk. As noted in lEspaillat et al. I (l2007bl . 1201 ll ) , the size of the 



hole is too large to be explained solely by the interaction of the binary with the disk. The 
outer regions of the circumbinary disk does not contri bute at the wavelengths observed (mid 



infrared), and for this reason, it is not included here (jD'Alessio et al. M2005I : lEspaillat et al. 



2007bl ). In principle, the emission of the wall is intrinsically variable with a period given by 
the orbital period of the binary system. This is because the two stars are changing their 
distance relative to different points on the wall. In addition, the wall has a different projected 
surface in the plane of the sky. However, in the case of CS Cha, the wall is l o cated too far 
from the stars, according to the SED modeling conducted by lEspaillat et al. I (1201 ll ). Using 
models incorporating two stars illuminating the wall, we found that the predicted difference 
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Fig. 1. — Comparison between the stellar spectrum used in lEspaillat et al. I ( 1201 ll ) and in 
this work. For reference the observations of CS Cha are represented wit h solid and pointed 
lines. The optical data i s shown as solid squares (jGauvin & Strom Ill992l ). The spectrum for 
the isolated star used in lEspaillat et al. I (120 111 ) is shown as a dashed line. The long dashed 
and point dashed lines represent the combined spectrum obtained from both stars for Case 
1 and 2, respectively. 



Table 1. Stellar parameters for the tested cases 



Case 


M +1 (M ) 






M^(M Q yr- 1 ) 


M* 2 (M ) 


n 2 (K) 




M, 2 {M Q yr- r ) 


1 


0.8 


4059 


1.74 


0.4£ - 8 


0.3 


3313 


1.29 


0.8E-8 


2 


0.9 


4204 


1.8 


OAE - 8 


0.1 


2998 


0.85 


0.8E-8 



en 

I 
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in the wall emission i s too small to explain the observations (1% at the 10/zm peak). As 



Espaillat et al. I (j201lf ). we propose that the mid-IR variability of CS Cha is due to variations 
in mass and illumination of optically thin dust in the circumbinary disk inner hole. We adopt 
a spatial distribution for this dust which it is based on results of numerical simulations in 
the literature. 



3.1. Outer wall 



We assume that the edge of the circumbinary disk is a vertical cylindrical wall illumi- 
nated by the stars. To the best of our knowledge, there are no calculations on the shape of 
the wall in circumbinary disks. For accretion disks around single stars, it has been proposed 
that the inn er wall is curved, beca use the sublimatio n temperature of the dust g rains depends 
on density (jlsella fc Natta 1120051 ) and on grain size (ITannirkulam et al. 1 120071 ) . However, in 
the case of a circumbinary disk, the inner disk truncation is produced by dynamical effects 
and not dust s ublimation. For the ca se of disks with embedd ed planets, there are detailed 
calculations by IVarniere et al. I (12004 ) and ICrida et al. I (120061 ) of the surface density profile 
of the gap, that in principle could be used to define a shape of the inner wall. However, the 
cases discussed in these papers correspond to low mass companions (planets) and it is not 
clear how one can scale this to binary systems. In addition, the density profile depends on 
viscosity, the disk aspect ratio, etc., which are unknown properties for CS Cha. 

We assume that the dust in each surface element of the inner edge of the circumbinary 
disk is in radiative equilibrium with the impinging radiation field from the stars. The flux 
received at each surface element of the wall is the addition of the contributions of both stars, 
taking into account the geometrical dilution of the stellar radiative flux, for the different 
distances of the stars at different orbital positions. Also, we assume that each star is accreting 
from a small circumste llar accretion disk, that receives mass from the circumbinary disk 
(jGimther fc Kiev 1120021 ) . It is usually assumed that for T Tauri stars surrounded by accretion 
disks, accretion shocks exist at the stellar surface that emit U V and optical radiatio n. We 
include the emission of these shocks as a source of irradiation (iMuzerolle et al. 1 120031 ). The 



luminosity liberated in these shocks is the accretion luminosity, ~ L acc , and we assume an 



effective temperature for the shocks of ~ 8000 K (ICalvet fc Gullbring I Il998l ). Thus, the 
spectral energy distribution of each star is a combination of their photospheric and shock 



emission. 



Th e zeroth and first mome nts of the radiative transfer equation are solved as in lCalvet et al. 



(Il99ll ) (ID'Alessio et al. II2005I ). The irradiation flux interacts with the dust in the wall, a 
fraction is absorbed and another fraction is scattered. The scattered fraction produces a 
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diffuse radiation field, which is finally absorbed at deeper layers. The transfer equation is 
integrated along the radial direction (parallel to the disk midplane) for each surface element 
of the disk inner wall and a temperature is calculated assuming radiative equilibrium. The 
temperature is a function of radial optical depth. In this calculation, we assume that the 
density of the wall at the edge of the circumbinary disk is high enough to allow thermal 
equilibrium between grains of different sizes and gas. Thus, a unique temperature describes 
the wall's thermal emission, for each wall location. 

For the calculation of the temperature on the surface of the wall [T(r = 0)], we use the 
optical properties of the dust at this temperature. The implicit equation for this temperature 
is then solved with an iterative procedure. The opacity evaluated at the wall surface temper- 
ature is used to calculate the temperature as a function of r. Finally, with the temperature 
radial distribution of each wall pixel, we calculate the emergent intensity, and taking into 
account the geometry of the wall and occultation effects, we calculate the contribution of 
the wall to the SED. 



3.2. Optically thin dust in the inner hole 



The models of gap formation in circumbinary disks show a decrease in mass surface 
density between an inner hole formed due to gravitational interactions of the binary and 
the disk (lArtymowicz fc Lubow II1994J ) and the outer disk, that suggests the gap is optically 
thin. Note that planets immersed in the disk ar e able to increase the size of the hole 
(iZhu et al. 1 1201 It iDodson-Robinson fc Salyk 1 1201 ll ) . Observations of ci r cumb inary disks at 



millimeter wavelengths, which resolve th e inner hole: iGuil 



Guilloteau et al. I ( 120081 ) for HH30, and I Andrews et al. I (120101 ) . con firm that it has a low 



oteau et al. I (Il999l ) for GG Tau 



density and l ow optical depth. Th is is the case for t ransitional disks : Hug 



for TW Hya, Hughes et al. I (120091 ) for GM Aur, and lAndrews et al. I (120111 ). The inner hole 



res et al. I (120071 ) 



might have optically thin d ust that contributes to the near and mid IR S EP, e .g., TW Hya, 
(Ukeson et al. Iboilh . and dCalvet et al~ll2005l : lispaillat et al. H2007bl Ja]. fioioh . 



We consider that the distribution of material inside the hole is not axisymmetric (see 
§ I3.2.ip . following the results from simulations of the interaction between the binary system 
and the disk. Also, since the dust in the hole is optically thin, with very low density, we 
calculate a different temperature for each grain size, assuming radiative equilibrium between 
the grain absorption of stellar and shocks radiation and emission of its own thermal radiation. 
Finally, we calculate the contribution to the SED of this optically thin material. 
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3.2.1. Dust spatial distribution 



The gravitational interaction between the disk and the binary system sets constraints 
on the possible configurati ons of the disk material. A binary system is able to cre ate a hole 
in the circumbinary disk (ILin fc Papaloizou I Il979l ; lArtymowicz &: Lubow I Il994l ). because 
there is an inner region without stationary orbital configurations. The stable orbits for the 
particles are located very close to each star (c ircumstellar disk s ) and further out, in tra- 



j ector ies around both stars (circumbinary disk; iPichardo et al. I (120051 ); iNagel fc Pichardo 



(120081 ) ). As one would expect, the closer the material is to one of the stars, the less im- 
portant its interaction with the other one. However, this configuration is perturbed when 
hydrodynamical interactions are taken into account. In this case, the material can lose an- 
gular momentum, and the resulting orbits end up connectin g their point of origin i n the 



circumbinary disk with each one of the circumstellar disks ( lArtymowicz &: Lubow 1 11996 



Giinther &: Kiev 1120021 ). In order to define a structure consistent with the results of these 
hydrodynamical simulations, we model the system as a circumbinary disk, two circumstellar 
disks and two streams of material, crossing the circumbinary disk hole, and connecting the 
outer circu mbinary disk with the inner cir cumstellar disks. Note that the larger hole present 
in CS Cha (jEspaillat et al. I l2007bl . 120111 ) cannot be explained only with the binary, thus, 
one requires another mechani sm to produce it , perhaps another stellar mass co mpanion or 
a multiple planetary system (jZhu et al. 11201 ll ; iDodson- Robinson fc Salyk 1120111 ). This will 
change the details of the structure in the hole, which we will not consider here. 



Artymowicz &: Lubow I (119961 ) present SPH simulations (IMonaghan 1119921 ) which follow 
the material that is traveling from the circumbinary disk to the stars. They note that 
the accretion is modulated in time with a periodicity given by the orbital period of the 
binary. The simple picture is that the material in the edge of the disk is perturbed near the 
apocenter of the binary, but requires a time delay to fall to one of the stars, arriving around 
pericenter. Of course the precise evolution de pends on the mass ratio and eccentricity of the 
binary system (lArtymowicz fc Lubow 1119961 ) . A result of this is that the time of evolution 
of the material through the gap is of the order of the orbital period (P=2482 days). This 
timescale is sma ll er than a typical dust c oagul ation timescale (~ 10 5 years for R ~ 10 AU, 
Weidenschilling I (ll977| ). iBirnstiel et al. I (120101 )). so we assume that the grain composition 
and size distribution do not change during travel. 

The circumbinary material in the outskirts of the disk initially evolves slowly, charac- 
terized by a viscous timescale evaluated at the radius of the circumbinary disk. Then, the 
material arrives to a radius defined by the gravitational intera ction of the binary system with 
the disk (lArtymowicz fc Lubow 111994 ; IPichardo et al. 1120051 ) . Starting from this radius, the 
matter lose its axisymmetry and flows in two streams. The two streams are launched from 
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two points called the saddle points. As noted in the last paragraph, the characteristic time 
for the evolution of the material in the gap is of the order of the binary orbital period. Thus, 
when the material arrives to the saddle points, it moves towards the inner circumstellar disks 
within a dynamical timescale. 

The saddle points are the analogous to the Lagrangian points in the 3-body circular 



restricted problem in classical mechanics ([Murray fc Dermott 



19991 ). For a circular binary, 



one can define a rotating coordinate system where the stars are at rest. Thus, there is 
a constant of motion, the Jacobi constant, and specifically two linear Lagrangian points 
( IMurray fc Dermott Ill999l ). that define the locus where material with enough value for this 
constant is able to move from orbits around the binary system to orbits connecting the 
circumbinary disk with circumstellar disks. In order to extend the analogy to an eccentric 
system, we note that the gravitational potential can be expanded as an infinite sum of terms. 
In many cases, one can describe the potential with only a few of them. For an eccentric binary 
system, the p otential can be simplified usin g the terms (m,/) = (0,0) and (m, /) = (2, 1) of 
its expansion (jArtymowicz fc Lubow 111 9941 ). Moreover, this potential rotates at a rate fi&/2, 
where fi& is the mean angular velocity of the binary. For this simplified case, two points 
can be defined (the corresponding Lagrangian points), which corotate with the potential. 
These are the launching points for the streams mentioned in the last paragraph, i.e. the 
saddle points. Analogous to the circular case, the material that accretes from the external 
regions of the disk is able to cross the radius where the saddle points are located, primarily 
at these locations. The saddle points are located at a radius R s = 1.587a, where a is the 
semimajor axis. Thus, we can conclude that at larger radii, the or bital backbone of the 
material are orbits around the binary system ( Pichardo et al. 12005). At smaller radii the 



material is located in tw o streams starting at the saddle points (jArtymowicz fc Lubow 111996 
Gimther fc Kiev 1 120021 1 and ending in the inner disks. 



The configuration used in this work, which is consistent with the orbital restrictions 
given by the binary system presence, is a ring of material between R s and R wa u, where R wa ii 
is the position of the inner edge of the outer disk. At R < R s , the material is located in the 
streams that connect the circumbinary and circumstellar disks. This is shown schematically 
in Figure [2J The ring corresponds to a region where the material evolves viscously, because 
this is consistent with the gravitational disk-stars interaction as described previously. The 
difference between this and the actual outer disk is that the ring is optically thin. As we 
mentioned above, in the case of CS Cha, for which we have Spitzer mid-IR fluxes, we do 
not include the contribution of t he circumbinary disk to the SEP, because it is too cold 
to emit in this wavelength range (jEspaillat et al. Il2007bl . l201ll ). Also, we assume that the 
circumstellar disks, streams and outer ring are all optically thin. 
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Fig. 2. — Sketch of the material configuration in the hole. It is not scaled. R s and R wa ii are 
the limits of the ring. The circumstellar disks and the streams are shown. 
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3.2.2. Dust temperature 



We assume the optically thin dust in the inner hole is in thermal equilibrium with the 
stellar radiation. However, the density in the gap is so low, that thermal equilibrium between 
gas and dust grains is not expected. Therefore, each grain with a different size would have 
a different temperature. 

The critical gas density for thermal equilibrium between dust grains an d gas can be esti- 

mated by assuming the collis ional timescale is equal to the thermal timescale ( 1 Chiang fc Goldreich 
19971 ; iGlassgold et al. 112004 ). The collisional timescale is given by 



1 



Tcol = Z n, W 

riHVira z 

where n# is the number density of the molecular hydrogen, v is the mean velocity of a 
molecule and a is the radius of the grain. We consider that the velocity is thermal, however, 
it is possible that there is a turbulent component, which amounts to a fraction of this velocity. 
Using the fac t that the typical velocity associa ted to turbulence in the q-parametrization o f 



the viscosity (IShakura fe Sunyaev et al. 1119731 ) in protoplanetary disks ( IHawley et al. Ill995[ ) 



is around 0.01 of the sound speed, we can conclude that including the turbulent velocity does 
not change our following conclusions. The thermal timescale can be estimated as the time 
that a particle of dust takes to radiate all its thermal energy at temperature T^(a), 



kTd(a) 



T the 



47raV 5s T rf (a)4' 

where k is the Boltzmann constant and <j,sb is the Stefan-Boltzmann constant. 



(2) 



Thus, the critical density is (IGlassgold et al. 1120041 ) 



Aa SB T d (af 
kv 



(3) 



and substituting typical numerical values, it can be written as 



(1.133 x 10 8 cm 



-3\ 



T 3 



1/2- 



(4) 



where T g is the gas temperature. 

Next, we would estimate the gas density associated with the optically thin material 
to compare it with this critical density. If we assume T g w T^, and take w 100-fT as 
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characteristic of the streams, then n cr = 3.34 x I0~ ll gcm~ 3 . Note that a change in T g of an 
order of magnitude corresponds to a factor of 3 change in n cr (see eq. [1]), which does not 
modify our conclusions. 

In order to get an order of magnitude for the density in the streams , we note that the area 
covered by the streams is around OAnR^ (lArtymowicz fc Lubow Ill996l ). We assume that the 
material falls from the inner edge of the ring towards the star in a few times the free fall time , 
which is 1.5 years. This is consistent with the fact that in lArtymowicz fc Lubow I ( 119961 ). 
the streams change on the order of an orbital period of the binary, T„ rh ~ 7 years. Thus, 
according to a mass accretion rate M = 10~ 8 M Q yr~ 1 (jEspaillat et al. Il2007bl ). the mass in 
the streams should be around 7 x 10~ 8 M^. Note that i n the hydrodynamical simulations of 



a disk around a binary system, iGiinther &: Kiev I (120021 ) assume a typical height in terms of 
the radius R in the disk, without dependence on the azimuthal angle, which is given by 



H(R) 



= 1,2 



ct\R-IL 







-1/2 



(5) 



where G is the gravitational constant, Mj is the mass of each of the stars, c s is the 
sound speed, and finally Ri is the distance of each star with respect to the origin. For Case 
1 and using a typical temperature of 100 K for molecular hydrogen and \R — Ri\ = R s , 
H(R) = 0.049.R S . Thus, a typical value for H is 0.1R S , which we consider here. Differences 
in H do not change the following conclusions. 

Given the mass in the streams, the area covered and the height calculated previously, 
the characteristic density in the streams is around n str = 10~ 14 gcm~ 3 , which means that 
n str << n cr . This implies that the gas and the dust grains are not in thermal equilibrium, 
and since the absorption coefficient of the grains depends on grain size, a, grains with different 
sizes have different temperatures, T(a). 

We divide the grain sizes in 100 bins equally spaced in log(a) between a min = 0.005/im 
and a max = 4/zm, distributed with a power-law profil e with exponent p = —3 .5. We take this 
value of a max according to the detailed modeling of lEspaillat et al. I (120111 ). For each size, 
we calculate the temperature from radiative equilibrium between the energy of the stellar 
and shock radiation absorbed by the grain and the energy lost by radiation. As expected, 
the temperature of the smallest grains is higher than the temperature of the larger grains 
(see discussion in § 15.21 and Figure [3]). We also estimate a mean temperature, as a unique 
temperature for the whole grain distribution, assuming the grains have a MRN power law size 
distribution, and use this temperature as a reference. This temperature is also calculated 
assuming radiative equilibrium with the radiation fields of the stars and their accretion 
shocks, but using mean absorption coefficients. As expected, at the same position in space, 
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the temperature of the smallest grains is higher and the temperature of the larger grains is 
lower than this mean temperature, due to a higher heating efficiency. Figure [3] illustrates 
this by showing temperature vs distance in the primary's stream (Case 1, GO observations), 
for the smallest grain (0.005/xm), the largest grain (4/im) and for a power law distribution 
of grain sizes between 0.005yum and 4/xm. 

Remember that the estimation of the temperature is done for all the locations where 
there is material and in each case, the temperature of each grain size is calculated. Th e 



grain composition adopted here is amorphous olivine, consistent with lEspaillat et al. 1 ( 120111 ). 
where 95% is olivine and 5% is crystalline silicates. In addition, we include organics and 
troilite. The abundances are ( s u = 0.0034, ( org = 0.001, and ( tro i = 0.000768, for silicates, 
organics and troilite, respectively. 



4. Parameters 

Since there are many input model parameters, we search for observational and theoret- 
ical constraints for some of them. In the next sections we describe the fixed and free model 
parameters. 



4.1. Fixed parameters 



For the best model for CS Cha, in lEspaillat et al. I (j2007b| ) and lEspaillat et al. I ( 1201 ll ). 
the inner edge of the circumbinary disk which we refer to as the "wall" is located at R wa ii = 
43 AU and R wa ii = 38AU, respectively. The diffe rence found in R v , a u. for CS Cha was 



due to using different dust opacities. The models of lEspaillat et al. I (120111) use a differen t 



distribution of silicates with respect to the models for CS Cha in lEspaillat et al. I ( )2007bl ). 
In the first modeling, 60% of the dust is silicates compared with 88% in the other case. For 
both cases, crystalline silicates correspond to a small amount. Here we use the most recent 
value, R W aii = 38AU, because it corresponds to a better fit of the 10/im band. At such a 
large distance, there is almost no difference between a model with one star and a model with 
two stars, because a « R wa ii (a = 3.703AU for Case 1, a = 3.587AU for Case 2). Thus, 



it is safe to assume the wall parameters to be the same inferred by lEspaillat et al. I ( 120111 ) , 
including the wall height (h = 7AU), and the inclination angle between the disk axis and 
the line of sight (i = 60°). 



As noted in § [21 we fix the value estimated by lEspaillat et al. I ( I2007bl ) for the mass 
accretion rate at 1.2 x lO _8 M ?/r _1 . We distributed this value: 0.4 x lO _8 M ?/r _1 and 
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Fig. 3. — Temperature along the stream connecting the primary star and the circumbinary 
edge, for the fit of the GO observation for Case 1. Small grains (solid line) correspond to 
a radius of a = 0.005/im, large grains (dashed line) correspond to a radius of a = 4/im and 
the grain size distribution (pointed line) is a power law of the form a -3 ' 5 , between the small 
and large grains 
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0.8 x 10 8 M ( r ) yr 1 , for t he pr imary and secondary, respectively. This is in agreement with 



Artymowicz fc Lubow I (119961 ). where the simulations with e = 0.1 show that the mass 
accretion rate to the secondary is a factor of two larger than the mass accretion rate to the 
primary star. Notice that the accretion luminosities are an order of magnitude lower than 
the star luminosity. As a result, we expect that the precise value of the former does not 
make a substantial difference in the resulting SED. The stellar parameters for Case 1 and 2 
are summarized in Table [TJ 



The streams start at the saddle points of the simplified potential: lArtymowicz fc Lubow 

( 119961 ) (see § 13.2. 1[) and end at each circumstellar disk. For simplicity their shape is a 
line that connects both positions. Although the thickness of the streams is a free param- 
eter, one condition should b e held: the flows fill around 10% of the surface of the hol e 
( lArtymowicz &: Lubow Ill996l ). Following the simulations of lArtymowicz &: Lubow I (119961 ). 
we can safely assume that the density of streams and outer ring are one tenth of the values 
associated with the compact circumstellar disks. We assume that both disks have a mini- 
mum radius R m in = 0.1AU and a maximum radius R max = 0.5 AU. As we mentioned above, 
the mass accretio n rate to the secondary i s abo ut two ti mes the value associated with the 
primary following lArtymowicz fc Lubow I (119961 ). see also iBate I (120001 ) . For the modeling, 
we translate this as a condition on the filling factor inside the hole associated to each stream. 
Thus, we assume that the area occupied by the primary star strea m is half the area of the 



stream reaching the secondary. The last requirement consistent with lArtymowicz fc Lubow 



(119961 ) is that the mass from the circumbinary disk that falls into any of the streams at 
apocenter is around twice the value at pericenter. Finally, in order to characterize the mass 
of the streams for every orbital location, we modulated the area of the streams (equivalently 
the mass) with a sinusoidal function in the azimuthal angle. The advantage of using this 
information to define the stream structure is that their geometrical parameters are fixed. 
The density on each stream is constant and it is fixed using the fit for the observed spectrum 
(see §E]). 



4.2. Free parameters 



Th e eccentricity is a param e ter th at is not restricted by the observations of lGuenther et al. 



(120071 ). lArtymowicz fc Lubow I (119961 ) point out that the variability of the streams mass de- 
pends on the eccentricity. The amount of material associated with the streams depends 
on the distance of the stars to the wall at apocenter, and this strongly depends on the ec- 
centricity. Note that for a circular binary, the minimum distance of each star to the disk 
does not change for the axisymmetric wall, thus the structure is stationary and there is no 
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variability in the amount of material in the hole. Thus, in order to explain the variability, 
the system must be eccentric. Also, the observed SED variatio ns suggest significant c hanges 



in the amount of material in the streams ( 17% according to lEspaillat et al. I ( 120111 )). For 
the modeling, in order to look for a best fit we take two values for the eccentricity, e = 0.1 
and e = 0.2. 

The remaining parameter needed to fully characterize the binary system is <fi. It is 
the angular position of the stars along the orbit at the time we are observing the system: 
= 0.57T corresponds to pericenter and = 1.57T to apocenter. The values we try in order 
to find the best fit for the earlier observation are 4>i — (0.0, 0.5, 1.0, 1.5)ir. The second epoch 
corresponds to 4>2, constrained by the time elapsed between the two observations, i.e., 1065 
days. We solve the equations of motion of the binary system, looking for configurations 
separated by this time. 

Finally, the dust mass surface density for the region assigned to the streams is a free 
parameter. The density is increased from zero until the flux level in the SED is reached. 



5. Results 

We calculate models at some representative points on the parameter space, as given in 
§ 14.21 For each pair of central stars (case 1 and case 2) and each pair of values (e, ^1,2), we 
calculate a SED adding the contributions of the stars, the wall and the optically thin dust, 
scaled by the dust mass surface density. The synthetic and observed SEDs are subtracted, 
and we calculate their corresponding x 2 ■ The mass surface density in each epoch is taken to 
be the one that minimizes \ 2 ■ 



5.1. Fitting the SED 

The fit is done independently for both epochs. In order to find the best fitting model, we 
search for the model that minimizes \ 2 ■ Because this is a small set of the parameter space, 
we cannot claim that we find a unique model. Due to the large number of assumptions 
and unknowns (for example, the exact configuration of the streams), degeneracy of the 
SED fitting is expected. Thus, one can only claim that we find models consistent with the 
observations. 

The fit for the GTO observation in 2005 is shown in Figure H] for Cases 1 and 2. The 
fit for the GO observation in 2008 is shown in Figure [5] also for both cases. The parameters 
for the best fit model for Cases 1 and 2 are the same, this is e = 0.1 and <fi = (0.5, 1.38)7r. 
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The shape of the streams' structure for the fit is shown in Figure [6] for Case 1, compared to 
the GTO and GO observation of CS Cha. The streams for Case 2 are almost identical to 
the ones of Case 1, thus, we do not show them here. 



5.2. Fitting the observed variability 

The fit in the last section is done for both epochs, but consistently taking into account 
the time difference between the observations. Thus, the variability can be explained with 
the previous fits. In Figure we show the Case 1 SED fits, for both observations, which 
are previously presented in Figures H] and Figured] is the same but for Case 2, showing 
the fits also presented in Figures H] and |5j We note that the best fit for either models with 
e = 0.1 or e = 0.2 is <p\ = 0.5. In order to check this tendency for every e, we ran models for 
e = 0.5 and e = 0.9. The tendency is corroborated, but the best models according to x 2 are 
worst for larger e. We ran models for e = 0.1 for other values of 0i, such that the whole set 
is 0i = (0.0,0.25,0.5,0.75, 1.0, 1.25, 1.5, 1.75). Even for this larger set, the best fit is found 
for (f) = 0.5. 

Given the observational cons traints and the proper ties we have fixed for the inner wall 



based on our previous modeling (jEspaillat et al. 1 120111 ) we find that any eccentricity that 



satisfies e > and <pi around 0.5, produces a configuration consistent with the observed 
SEDs in the two epochs. This means that the eccentricity e is not constrained by our models 
and it appears that only observations of the detailed orbit of the binary system would allow 
to quantify this parameter. 

The SED of the present mod els (for both, Ca s e 1 an d Case 2) are similar to the SED of 



a single star model presented by lEspaillat et al. I ( 1201 ll ). In the single star model the SED 



variation is related to a change in the mass of the optically thin dust in the inner hole. In 
the present model with two central stars, we also consider a change in dust mass, but besides 
this there is a source of variability given by the change in illumination produced by the stars, 
which is absent in a single star model. 

As we show, the variability can be explained by the differences for each orbital config- 
uration of the binary system, in the illumination of the material in the circumstellar disks 
and streams. Also it is important the amount and the geometrical configuration of the dust 
in terms of the orbital configuration. A visualization of this effect is given in Figure |H] where 
for the Case 1 fit, the intensity for a wavelength of 10/zm is presented. One can clearly see 
that the contribution to the emission is larger for the configuration associated to the GO 
observation (Figure E]). Notice that in Figure |H1 the variation in area for the streams means 
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X(jum) \(/zm) 

Fig. 4. — The best fit for CS Cha for Cases 1 (left panel) and 2 (right panel) (see Table [1]) for 
the GTO observation on 2005. The solid line is the SED observation. The pointed line is the 
hole SED for the model. The dashed line corresponds to the wall flux. The addition of both 
stellar SEDs are represented with a point-long-dashed line. Finally, the small-long-dashed 
line is the total flux for the models. 
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X(/zm) \(jum) 



Fig. 5. — The best fit for CS Cha for Cases 1 (left panel) and 2 (right panel) (see Tabled]) 
for the GO observation on 2008. The lines have the same meaning as Figure HI 
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Fig. 6. — Diagram of the structure of the streams viewed pole on, for the best fit for the Case 
1 (see Tabled]) for the GTO observation (upper plot) and the GO observation (lower plot). 
The observer is at an inclination of 60°. The crosses (+) represent the primary star and the 
rotated crosses (x) represent the secondary star. The region inside both lines corresponds 
to the material associated to the streams. The circumstellar disks are drawn and the label 
'X' refers to the saddle points, where the streams are launched. The surrounding circle 
represents the inner boundary of the optically thin ring of material. 
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0.9 1 1.1 1.2 1.3 

log \(/im) 

Fig. 7. — The modeled variability for CS Cha for the Case 1 (see Table [1]). The solid lines 
represent the SED for the GTO observation (2005) and the GO observation (2008). The 
dashed and the pointed lines are the modeled SED for both epochs. These correspond to 
the models shown in Figures H] and [H 



-23 - 




0.9 1 1.1 1.2 1.3 

log \(/im) 

Fig. 8. — The modeled variability for CS Cha for the Case 2 (see Table [1]). The solid lines 
represent the SED for the GTO observation (2005) and the GO observation (2008). The 
dashed and the pointed lines are the modeled SED for both epochs. These correspond to 
the models shown in Figures H] and [H 



-24- 



a larger amount of mass, and the gray scale corresponds to changes in the illumination. 



The associated dust mass in circumstellar disks and streams is around 10 12 M Q , of the 



same order of the estimate by lEspaiUat et al. I (1201 if ), where the material was located only 
in one circumstellar disk. The dust surface density on each stream is uniform and can be 
calculated with the total mass of dust and the area associated to each stream. In Table [3], we 
present the dust surface density for the primary and the secondary streams, for Case 1 and 2, 
and for the configurations associated with both epochs. In order to estimate from the model 
a mass accretion rate we do the following. We assume that the gas attache d to the dust in 



the disk s and streams is responsible for the mass accretion rate estimated in lEspaillat et al. 
(l2007bl ). Using the dust abundances, the mass of gas is of the order of 1O~ 1O M , if we now 



use a free fall time from the saddle point as the fastest time for the accretion, the mass 
accretion rate i s around 10~ 10 M^yr -1 . This value is two orders of magnitude smaller than 



the estimate in lEspaillat et al. I ( 12007bl ). A possible explanation for this is that the dust to 
gas ratio is lower than usually expected . A way to do th is is through a filtering effect in 
the inner edge of the circumbinary disk. iRice et al. I (120061 ) note that the pressure gradient 
in this location act as a filter, only allowing small dust grains to trespass it, holding the 
larger ones. However, the typical critical size is around 10/im, which is larger than the 
maximum value taken in our work, a max = 4/xm. Another idea is that the mass accretion 
rate inferred from U-band excess has uncertainties due to the estimation process in itself 
and due to the fact that the system is a binary. An improvement of this estimate requires a 
detailed description of the evolution of the particles (solid and gaseous) si nce they leave the 



inner edge of the circumbinary disk to finally accrete at one of the stars (jGimther &: Kiev 



20021 ). Noteworthy, an U-band excess is not always an indicative of accretion, because one 
has to check the sha pe of the lines profiles to safely assume that the system is accreting 
(jlngleby et al. 1120111 ). However, a complete study able to elucidate the dynamics of the hole 
material is beyond the scope of this paper. 



In previous works ( ID'Alessio et al. Il2005t iNagel et al. Il2010l ; lEspaillat et al. 11201 ll ). the 
dust temperature in the hole and the emissivity are grain-size-mean values calculated using 
an opacity, given by the average opacity over the grain size distribution. In here, since the 
dust in the hole has suc h a low density and cannot be considered in ther mal equilibrium 
with the gas (see § 13.2.21 (jChiang fc Goldreich Ill997l ; iGlassgold et al. 1120041 ). we calculate a 
dust temperature and emissivity as a function of grain radii. For the sake of evaluating the 
differences between both approaches, we calculate a model with the grain-size-mean opacity 
to fit the SED. We find that the resulting SEDs are very similar, but for the model with the 
grain-size-mean opacity assumption, the mass of the optically thin material is 19% of the 
mass of the models where no thermal equilibrium is assumed. This difference between both 
approaches means that it is worthed to do the proper estimate. 
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Intensity map 




-1.5 -1 -0.5 0.5 1 1.5 

x 

Fig. 9. — Intensity map at A = 10/zm for the Case 1 fit. The gray scale shows the log of the 
intensity, which is presented at the right with units erg cm~ 2 s~ 1 ster~ 1 . The coordinates x 
and y are given in units of a. We made a superposition of both configuration (see Figure [6]). 



Table 2. Orbital parameters for the cases tested 



Case n( ATI) 




< r\ n ) 


1 3.703 


0.1 


0.0 






0.5 






1.0 






1.5 




0.2 


0.0 






0.5 






1.0 






1.5 


2 3.587 


0.1 


0.0 






0.5 






1.0 






1.5 




0.2 


0.0 






0.5 






1.0 






1.5 



Table 3. Surface density (E) of the streams 



Case 


Epoch 


Stream to 


E(gr cm 2 ) 


i 


V i 1 V / 


jji iiiidi y 






GTO 


secondary 


1.92E-5 




GO 


primary 


1.58E-5 




GO 


secondary 


1.45E-5 


2 


GTO 


primary 


1.91E-5 




GTO 


secondary 


4.05E-5 




GO 


primary 


2.85E-5 




GO 


secondary 


1.43E-5 
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6. Summary and conclusions 



We have shown that the Spitzer IRS spectrum of the binary CS Cha is consistent 
with the emission from the inner disk structur e generated by the do uble system. This is 



a step forward with respect to the modeling in lEspaillat et al. I ( 1201 ll ). because we include 



the binary system and more details about the hole dust configuration. The presence of a 
binary means that the disk is gravitationally truncated, forming an inner edge (wall) which 
is directly illuminated by the stars. This would be a good candidate to follow up for mm 
imaging to confirm the hole size. 



As already discussed by lEspaillat et al. I (120111 ). the mid-IR variability of CS Cha is 
restricted to the 10/im silicate band. This implies that it is not produced by the variable 
illumination of the wall, which has a large radius as inferred from SED fitting, but by 
variations in the emission of optically thin dust in the inner hole of the circumbinary disk. 
We adopt a mass distribution inside the hole consistent with dynamical simulations (see next 
paragraph), i.e., the optically thin dust is located in a ring, two streams and two circumstellar 
disks. The mass of the streams is modulated by the orbital configuration of the stars, which 
combined with a variable illumination from both stars, can explain the observed variability 
for a reasonable set of parameters. 

The distribution of the material in the hole is consistent with the dynamical restric- 



tions due to the presence of a binary system (lArtymowicz fc Lubow I Il994j ; iPichardo et al 



20051 ). The material is located in a ring, streams and circumstellar disks (see § I3.2.ip . This 
configuration is able to explain the observed variability. A SED fit is done for the optically 
thick regime using hydrodynamical simulations of th e form ation and evolution of the hol e 
material in a circumbinar y disk in iGiinther &: Kiev I (120021 ) and iGimther &: Kiev I ( 12004 ) . 
de Val-Borro et al. I (120111 ) do similar simulations for the close binaries V4046 Sgr and DQ 
Tau, aiming to explain the shape variability of gas lines. In this case, the density of the 
structure points to the optically thick regime. Thus, it is important to point out that for the 
first time (as far as we know), the SED fitting of the mid-infrared variability, for the emis- 
sion of optically thin material inside the hole of a circumbinary disk, using a configuration 
consistent with the theory and simulations is successfully done. 

The masses for the esti mated optically thin du st are of the same order of magnitude 
than the values estimated in lEspaillat et al. I (1201 if ), however, note that the configurations 
of the material in the gap are quite different between both cases. Given the unknowns, we 
decide to calculate a small set of models in order to illustrate that they can explain the 
observed SEDs and variability, but we do not claim this models are unique. However, new 
observations of this system will be able to increase the constraints, for a more extensively 
testing of our model. In the future, this model will be tested for more binary stars with 
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multiple observations in the infrared. Unfortunately, by now the cryogenic mission on board 
the Spitzer Space Telescope is over. The ongoing warm mission is just able to confidently see 
at the smaller range of wavelengths. Thus, we require the launching of the James Webb Space 
Telescope to pursue this project further. The interferometer ALMA can see at a wavelength 
(330/xm) which is at least an order of magni tude larger than the Spitzer wavelength window. 
At this frequency, I Wolf fc D'Angelo I ( 120051 ) study simulations images of a planet embedded 
in a protoplanetary disk. They conclude that the planet with an accretion luminosity of 
10~ 5 L* is able to perturb the disk in order to be detected by ALMA. In our case, the ratio 
of the emission from the stellar peak and the estimated emission at 330/im is around 10 6 , 
thus, we think that the streams can barely be seen with ALMA. 
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